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ABSTRACT 


The  method  of  two  loci  for  scallng-up  the  drag  reduction 
by  additives  from  measurements  In  not-too-small  diameter 
pipes  to  large  diameter  pipes  Is  extended  to  small 
diameter  pipes.  This  Is  accomplished  by  Including  the 
effects  of  the  viscous  sublayer  and  buffer  layer  In  for¬ 
mulating  a  prediction  chart.  A  mixing-length  model  proved 
more  accurate  than  three  other  models  of  the  viscous 
sublayer  and  buffer  layer  In  a  comparison  of  predictions 
with  measured  data.  The  prediction  chart  may  also  be  used 
to  obtain  a  slmllarlty-law  drag-reduction  characterization. 


INTRODUCTION 

The  prediction  of  drag  reduction  by  additives*  In  large  diameter  pipes  from 
measurements  In  small  diameter  pipes  using  the  method  of  two  loci^  has  been 
shown  by  Sellln  and  Ollls^  to  be  accurate  when  the  diameters  of  the  test  pipes 
are  sufficiently  large.  From  an  engineering  viewpoint.  It  may  be  more  economi¬ 
cal  to  use  small  pipes  to  measure  the  actual  drag  reduction  characteristics  of 
commercially  available  additives  which  may  vary  In  molecular  sizes  due  to  the 
manufacturing  process.  It  Is  the  Intention  of  this  paper  to  extend  the  method 
of  two  loci  to  small  pipes  by  Including  the  effects  of  the  viscous  sublayer  and 
the  buffer  layer  which  %fere  not  Included  In  the  original  method  developed  for 
large  pipes  ^ . 

The  problem  of  scallng-up  the  results  from  smaller  pipes  to  larger  pipes 
arises  from  a  two-parameter  variation  of  the  friction  factor  with  not  only  the 
usual  Reynolds-number  dependence  but  also  a  nondlmenslonal  diameter  for  the 
same  concentration  and  species  of  additive^.  Fortunately,  the  classical  velo- 

3  If 

city  similarity  laws  ’  provide  a  characterization  of  drag  reduction  which  Is 


*  Host  of  the  material  In  this  paper  was  presented  to  the  Third  International 
Conference  on  Drag  Reduction,  July  2-5,  1984  at  the  University  of  Bristol, 
England . 


independent  of  pipe  dlaaeter.  The  original  nethod  of  two  loci  considered  the 
logarithmic  similarity  law  to  be  valid  up  to  the  wall,  thus  Ignoring  the  effects 
of  the  viscous  sublayer  and  buffer  layer.  The  method  was  hence  applicable  to 
large  pipes  only.  By  Including  the  effects  of  the  viscous  sublayer  and  buffer 
layer,  the  nethod  of  two  loci  can  be  extended  to  small  pipes.  Four  existing 
models  of  the  similarity  laws  for  the  viscous  (laminar)  sublayer  and  the  buffer 
layer  are  to  be  considered,  namely, 

1.  A  variable  thick  laminar  sublayer  and  an  associated  buffer  layer^ ; 

2.  A  constant  laminar  sublayer  and  an  elastic  layer  In  lieu  of  the 
buffer  layer’; 

^  7 

3.  A  mixing-length  formulation  *  ;  and 

4.  An  extended  viscous  sublayer  . 

A  brief  review  of  the  velocity  similarity  laws  Is  presented.  The  four 
models  are  used  to  obtain  four  friction-factor  formulations.  The  results  are 
presented  graphically  fo"  the  extended  nethod  of  two  loci.  Predictions  from  the 
four  models  are  compared  with  experimental  data.  An  appendix  shows  how  to 
determine  a  simllarlty-law  drag  characterization  from  small-diameter  pipe  flows. 

VELOCITY  SIMILARITY  LAWS 

In  general,  the  inner  similarity  laws  for  drag-reducing  additives  (3)  may 
be  stated  nondlmenslonally  as 


UL 


Uxl 


(1) 


where  u  >■  u/u^;  u  ••  mean  turbulent  velocity  parallel  to  the  pipe  wall;  u^ 


✓t^/p;  *  wall  shearing  stress;  p  »  density  of  fluid;  y  *  u^y/v;  y  *  nor¬ 


mal  distance  from  the  wall;  v  ~  kinematic  viscosity  of  fluid;  i  *  characteristic 

A 

length  representing  the  additive;  C  •  concentration  of  additive;  A  >  species  or 

A 

kind  of  additive.  The  characteristic  length  Z  In  this  paper  serves  mostly  as  a 
dummy  variable. 

A  characteristic  time  t  or  mass  m  may  be  substituted  for  the  characteristic 
length  t  where  i  -  /tv  *  (mp)^^^. 

The  outer  similarity  law  which  does  not  directly  depend  on  the  presence  of 
additives  may  be  stated  for  fully  developed  pipe  flow  as 


where  U  ■  velocity  at  the  center  of  the  pipe  and  r  *  radius  of  the  pipe. 

The  overlapping  of  the  Inner  and  outer  laws  results  In  logarithmic  laws 


u*  “  A  in  y*  + 


where 


6.= 


may  be  termed  the  drag-reduction  characterization.  Another  drag-reduction 


characterization  Is 


A6  =6.-%,,  =/ 


where  q  =  value  of  In  the  absence  of  the  additive.  A  and  ^2  constants 
for  pipe  flow. 

If  the  logarithmic  velocity  law.  Equation  (3)  is  considered  to  extend  to 
the  pipe  center,  a  correction  has  to  be  applied.  The  nondlmensional  velocity  u* 
across  the  pipe  Is  to  be  represented  by  two  components 


*  * 

u  “  +  U2 


where  u^*  Is  the  contribution  of  the  laminar  (viscous)  sublayer,  the  buffer 
layer  and  the  logarithmic  layer  extended  to  the  pipe  center. 


Uo  is  the  correction  represented  by 


where  w  is  the  Coles  wake  function  and  q  Is  the  wake  modification  function^. 
2 


Polynomial  fits  to  w/2  and  q  are  given  by 


(9) 


and 


1  - 


*  x3 


(10) 


where  r  =  u^r/v 

PIPE  FLOW 

For  the  fully-developed  flow  of  drag-reduction  solutions  through  pipes  the 
Fanning  friction  factor  f  may  be  stated  functionally  as 


f  =/|Ir,  D/1,  C,  a  ] 


(11) 


2 

where  f  =  Fanning  friction  factor  =  2t^/pV  ;  V  =  average  velocity  across  the 

pipe;  R  =  pipe  Reynolds  number  =  VD/v;  D  =  diameter  of  the  pipe  =  2r. 

A 

For  the  same  concentration  C  and  species  A  of  additive,  the  friction  factor 
varies  not  only  with  Reynolds  number  R  but  also  with  nondlmenslonal  diameter 
D/l;  hence  the  scale-up  problem.  Another  commonly  used  friction  factor  Is  the 
Darcy-Welsbach  friction  factor  \  where 


X  =  4f 


(12) 


The  average  velocity  V  across  a  circular  pipe  may  be  determined  from  the 
simllarlty-law  velocity  profile  by  means  of  ^ 

^  ■  -f.  (/“■  -  r.  ly  >■ 

which  is  derived  in  Reference  3. 

Also  by  definition 


(14) 


♦Note  that  hydraulic  engineers  seem  to  prefer  the  Darcy-Welsbach  friction  factor 
using  either  the  symbol  f  or  X.  Chemical  engineers  seem  to  prefer  the  Fanning 
friction  factor  using  the  symbol  f  as  used  here. 


and 


which  bridge  sinllarlty-law  variables  and  pipe  variables. 

EXTENDED  METHOD  OF  TWO  LOCI 
If  the  logarithmic  velocity  law.  Equation  (3),  is  assumed  to  hold  up  to  the 
wall,  that  Is,  without  any  effect  from  the  viscous  sublayer  and  the  buffer 
layer,  then  by  means  of  Equations  (7),  (8),  (9),  and  (10)  there  results 


If 


(16) 


For  ordinary  fluids  (AB  ~  0),  this  equation  has  the  same  form  as  the  well- 
known  Prandtl-Karman  formula  for  pipe  flow 


—  =  4-.0  ^  -  0. 4- 


(17) 


Now  for  the  same  concentration  C  and  species  A  of  additive  the  similarity- 
law  characterization  AB  is 


AB 


A  n 


(18) 


Hence  the  values  of  both  AB  and  have  to  be  satlslfled  for  pipes  of 


different  diameters. 


It  is  obvious  that,  for  a  plot  of  1/  /f  against  /T  R,  lines  of  constant  AB 
are  parallel  to  each  other  in  accordance  with  Equation  (16).  This  Is  the  first 
locus  which  only  applies  to  sufficiently  large  pipes  where  the  effects  of  the 
viscous  sublayer  and  the  buffer  layer  are  negligible. 


The  second  locus  has  to  satisfy  the  saae  value  of  u^l/g,  idilch  by  definition 


For  the  sane  concentration  and  species,  1  Is  constant.  Thus  In  going  from  pipe 

A 

diameter  to  pipe  diameter  D2,  for  the  same  value  of  u^l/v. 


.(R  1%) 


This  represents  the  second  locus  which  does  not  exclude  the  effects  of  the 
viscous  sublayer  and  the  buffer  layer. 

For  smaller  pipes  where  the  effects  of  the  viscous  (laminar)  sublayer  and 
the  buffer  layer  are  not  negligible,  the  lines  of  constant  &B  on  a  l//f  -  ZfR 
plot  are  no  longer  straight  and  parallel.  Hence  the  curved  lines  of  constant  AB 
have  to  be  plotted  at  sufficiently  close  Intervals  so  that  the  locus  of  constant 
AB  Is  drawn  from  a  measured  point  by  Interpolating  between  two  nearby  lines  of 
constant  AB.  This  is  shown  In  Figure  1. 

Four  analytical  models  of  the  viscous  (laminar)  sublayer  and  the  buffer 
layer  are  now  to  be  examined  in  order  to  determine  the  curved  lines  of  constant 


THICK  LAMINAR  SUBLAYER 

The  first  analytical  model  proposed  for  the  buffer  layer  seems  to  be  that 
of  reference  3  which  was  based  on  an  eddy  viscosity  model.  The  result  Is  a 
laminar  sublayer  which  thickens  with  increasing  values  of  AB. 

As  usual  the  laminar  sublayer  of  thickness  y^*  Is  given  by 

-  y  .  0  <  y  <  y^  (21) 

Outside  the  laminar  sublayer  a  logarithmic  relation  holds  for  both  the  buffer 
layer  and  the  logarithmic  layer  which  Is 

Uj^*  -  A  tn(y*  -J)  +  Bj^,  y*  >  y^^*  (22) 

This  relation  joins  the  laminar  sublayer  In  both  value  and  derivative.  The 
result  Is 

J-Bj-A+AinA  (23) 


+  A  In  A 


6 


and  for  ordinary  fluida 

Jt  -  Bi  0  -  A  -f  A  in  A  (25) 

Deterainlng  the  average  velocity  from  Equation  (13)  and  converting  to  /f 
and  /Ir  coordinates  yields 

=  (26) 

where 

4,  «  2,  KIJo  ^ ♦  2(Xo-*  Aft^  ♦  A  -^  tS^ 

It  is  to  be  noted  that  this  equation  is  the  same  as  that  of  Equation  (16) 
except  for  the  added  terms  due  to  the  laminar  sublayer  and  the  buffer  layer 
which  decrease  as  /fR  Increases  for  larger  diameter  pipes. 


ELASTIC  LAYER 

5 

On  the  basis  of  a  line  of  maximum  drag-reduction  in  pipe  flow.  Virk  et  al 
proposed  a  logarithmic  relation  for  the  buffer  layer,  which  was  renaited  the 
elastic  layer,  of  the  form 


*  'T'  ,  *  'Sr' 

u  »  A  in  y  -  B  (27) 

where  A  and  B  are  constants.  This  relation  begins  at  the  intersection  of  the 

A 

logarithmic  line  and  laminar  line  for  the  ordinary  fluids,  y^^  .  As  stated,  no 
buffer  layer  is  allowed  for  ordinary  fluids  so  that  y-j^  is  a  constant. 

For  drag-reducing  fluids  the  elastic  line  ends  at  y*^^  which  is  the  intersec¬ 
tion  of  the  elastic  line  and  the  logarithmic  line  or 


%  * 


^*<o  f  0  •»’  A  & 


(28) 


Performing  the  required  integrations  of  Equation  (13)  results  in 


IT 


.A _ 

IT  ffa  Ui&f 


(29) 


where 


j .  t  -  i  A  'i*  j 


And  m  3  *** 

Recently  this  elastic-layer  fomulation  was  used  in  an  iterative  calculation 
method  by  Matthys  and  Sabersky^^  in  scaling-up  pipe  data. 


MIXING  LENGTH 

Poreh  and  Dimant^*^  extended  the  Prandtl-van  Driest  mixing-length  formula¬ 
tion  to  drag-reducing  solutions  as  follows.  Both  the  viscous  sublayer  and 
buffer  layer  are  included. 

According  to  mixing-length  theory  the  turbulent  shear  stress  is  related 
to  the  velocity  gradient  as  2, 

»  C  I  — \  (SO) 


where  X  is  the  mixing  length. 

Including  the  laminar  shear  stress  results  in  nondimens ional  terms  as 


where  t*  -  t/t^,  1*  ■  u^l/u  and  t  ■  shear  stress  in  fluid. 

Solving  as  a  quadratic  quation  yields  a  velocity  profile 

The  Prandtl-van  Driest  mixing-length  fonmilatlon  is  nondlmenslonally 


where  k  ~  von  Karman  constant  «  1/A;  \  >  u^X/u  and  X  *  van  Driest  factor. 

For  pipe  flow  the  shear  stress  is  given  as 

*  ,  * 

T  ■  1  -  ^ 


8 


The  van  Driest  factor  X  nay  be  correlated  with  the  drag-reduction  charac¬ 
terisation  or  &B  as  follows. 

At  high  values  of  r*  where  r*  "  1  and  where  the  wlxlng-length  velocity  pro¬ 
file,  Equation  (32),  nerges  Into  the  logarithmic  velocity  profile 


e.  = 


i/  -  K  Ik-V" 


A  A 

where  y  Is  a  sufficiently  large  value  of  y  .  B^  Is  thus  a  function  of  X  . 

If  Equation  (32)  Is  considered  to  extend  to  the  pipe  center,  u*  becomes  Uj^* 
and  the  outer  correction  U2*,  as  assumed  by  Poreh  and  DlsMnt,  Is 


The  evaluation  of  the  average  velocity  bv  Equation  (13)  has  to  proceed 
numerically  In  order  to  develop  a  1//F  vs.  /Tr  relationship.  Tlederman  and 
Relschman^ ^  show  that  the  Prandtl-van  Driest  mixing-length  formulation 
Indirectly  fits  the  measured  velocity  profiles.  The  eddy  viscosity  relationship 
actually  used  by  them  may  be  shown  to  be  based  on  the  Prandtl-van  Driest  mixing 
length. 

EXTENDED  VISCOUS  SUBLAYER 

An  extended  viscous  sublayer  proposed  by  Nasan  for  ordinary  fluids  was 

D 

applied  to  drag-reducing  fluids  by  Kwack  and  Hartnett  .  The  effect  of  an  eddy 
viscosity  was  added  to  the  laminar  flow  to  give  a  velocity  profile  for  both  the 
viscous  sublayer  and  the  buffer  layer  such  that 

A  A  A**  A® 

u  »y  +Cj^y  +C2y  (37) 

where  c^  and  C2  are  constants  to  be  determined.  The  fourth  power  of  y*  is  chosen 
to  provide  a  third-power  variation  of  eddy  viscosity  with  y  at  the  wall.  The 
Junction  of  Equation  (37)  with  the  logarithmic  law.  Equation  (3),  is  at  y^. 

A 

The  values  of  Cj^,  C2  and  y^  vary  wth  Bj^  or  &B  and  are  determined  by  con¬ 
tinuity  at  the  junction  y  so  that 
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Using  Eq.  (13)  and  converting  to  and  /7r  coordinates  results  In 


jr  10 M*?!'  jfti 


where 


and 
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NUMERICAL  AND  GRAPHICAL  RESULTS 


Charts  for  Method  of  Two  Loci 

The  charts  for  the  extended  nethod  of  two  loci  are  shown  In  Figures  2,3,4, 
and  S  for  the  thick  lanlnar-aublayer  nodel,  the  elastic-layer  nodel,  the  nlxlng- 
length  nodel,  and  the  extended  viscous-sublayer  model  respectively. 

For  the  case  of  ordinary  fluids,  AB  »  0,  the  usual  Prandtl-Karman  fomula. 
Equation  (17),  Is  plotted.  Values  of  A  -  2,5,  •  5.5,  A^-  11.7,  17.0 

are  used  In  the  calculations  for  the  thick  laminar-sublayer  nodel,  the  elastic- 
layer  nodel  and  the  extended  viscous  sublayer  nodel.  A  value  of  »  11.6  Is 
used  for  the  elastic-layer  nodel. 

The  chart  for  the  nlxlng-length  model  is  developed  from  a  f;raph  of  computed 
results  presented  by  Dimant  and  Poreh  (7). 

Coiaparlson  of  Predictions 

The  experinental  data  of  Sellln  and  Ollis  (2)  for  the  drag-reduction  in  a 
5  nn  diameter  pipe  and  in  a  50  nn  diameter  pipe  are  used  to  compare  predictions 
from  the  three  methods  in  scaling-up  the  data  from  the  Smn<  to  the  50  mm  pipe. 

The  measured  data  are  for  a  10  wpm  concentration  of  4  I'fimmercial  polyacrylamide, 
Alcomer-llOL.  Solid  symbols  represent  experiments'  dsis  while  )pen  symbols 
represent  predictions. 


As  shown  In  Figures  2,  3,  4  and  5»  the  nlxlng-length  wodel  gives  excellent 
predictions  for  the  few  data  presented.  The  thick  laalnar-sublayer,  the  elastic- 
layer  and  extended  viscous-sublayer  models  show  a  poor  comparison  between  pre¬ 
dicted  and  measured  data.  This  seems  to  Indicate  that  the  mixing-length  model 
gives  the  best  representation  of  the  viscous  sublayer  and  buffer  layer. 

CONCLUDING  REMARKS 

The  Inclusion  of  the  effects  of  the  viscous  (laminar)  sublayer  and  the 
buffer  layer  definitely  seems  to  Improve  the  accuracy  of  the  original  method  of 
two  loci  especially  for  small-dlameter  pipes.  This  may  be  seen  by  noting  the 
curved  lines  of  constant  AB  In  the  l//f  vs.  /fR  plot  which  deviate  from  the 
straight  lines  specified  In  the  original  method  of  two  loci.  At  higher  values 
of  td\lch  apply  to  large  diameter  pipes,  the  lines  of  constant  AB  start  to 

become  straight  and  parallel  In  accordance  with  the  original  method  of  two  loci. 
This  explains  the  success  of  Sellln  and  0111s  (2)  In  making  predictions  from  25 
to  50  mm  diameter  pipes  and  the  lack  of  success  In  making  predictions  from 
smaller  than  25  mm  diameter  pipes. 

The  van  Driest  mixing-length  model  gives  predictions  In  scallng-up  the  data 
from  5  mm  diameter  pipe  to  a  50  mm  pipe  which  agree  well  with  the  available 
experimental  data.  Hence  the  van  Driest  model  as  given  In  Figure  4  Is  recom¬ 
mended  for  the  extended  method  of  two  loci. 


^  m  »  •  •  m  • 
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APPENDIX 


DETERMINATION  OF  SIMILARITY-LAW  DRAG-REDUCTION  CHARACTERIZATION 


The  AB-characterizatlon  may  be  also  determined  from  the  charts  l/7f  vs.  /fR 

A 

and  AB^  as  a  function  of  u^l/v  for  the  particular  concentration  and  type  of  addi¬ 
tive  Involved. 

The  test  point  (l//f.  /fR)  for  a  particular  diameter  pipe  D  gives  the  AB 
value  from  an  Interpolation  between  the  lines  of  constant  AB. 

From  definition,  the  value  of  u^f  Is  given  here  as 


u^l  -  ^  /fR  .  ^  (Al) 

“D"  2  D 


A 

In  many  applications,  the  characteristic  length  I  serves  as  a  dummy  quantity. 

AB  characterlz:.tlon  may  then  be  used  to  determine  the  viscous  drag  of  bodies  by 
a  boundary-layer  calculation. 
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FIGUHE  4  -  Method 


Method  of  Two  Loci  for  Extended  Viscous-Sublayer  Model 


DTNSROC  ISSUES  THREE  TVPES  OF  REPORTS 


1.  DTNSROC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY.  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


